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ABSTRACT 

The measurements of the possible gravitational microlensing events are analysed with 
a simple yet accurate disc-halo model of the Milky Way Galaxy. This comprises a 
luminous exponential disc embedded in a flattened dark matter halo with density 
varying like distance" ' . Including a disc has the important effect of lowering the 
implied masses of the dark matter objects. For the possible detection reported by 
Alcock et al (1993), the inferred mass of the lens lies in the range ^ 0.01 — O.15M0. 
The candidate events of Aubourg et al (1993) have slightly larger implied masses of 
~ 0.025 — 0.35 Mq and ~ 0.03 — 0.45 Mq respectively. These are consistent with 
the deflecters being either brown dwarfs or low mass stars. If there is no disc dark 
matter and the halo is completely composed of baryonic dark objects of typical mass 
^ 0.08 Mq, then the monitoring of 1.8 x 10^ stars in the Large Magellanic Cloud will 
provide at least 8 — 9 detections a year, in the limit of 100% efficiency. 
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1 INTRODUCTION 

Paczyriski (1986) made the influential proposal that dark ob- 
jects - possibly brown dwarfs, low mass stars or black holes 
- in the halo of the Milky Way can act as gravitational mi- 
crolenses, causing the occasional amplification of the images 
of stars in the Magellanic Clouds. He used the simplest 
possible model of a galactic halo - the singular isothermal 
sphere - to estimate the order-of-magnitude efi'ects. Griest 
(1991) refined the model by introducing a core radius and 
analysing the importance of relative motion of observer and 
source. Two groups (Alcock et al 1993; Aubourg et al 1993) 
have now claimed possible detections of microlensing events 
of Large Magellanic Cloud stars. Evidently, this warrants 
the development of a more accurate representation of the 
Milky Way galaxy to enable a sophisticated confrontation 
with the observations. Such is the purpose of this Letter. 



2 A DISC- 
WAY 



HALO MODEL FOR THE MILKY 



The observables that can be measured by the experimen- 
tal programs are the microlensing rate F and the average 
duration of events (fe). To compare the predictions of a 
theoretical model with the data is only possible if the veloc- 
ity distribution of the massive compact halo objects (hence- 
forth Machos) is known. As the Machos are coUisionless, 
the distribution function must obey Jeans' theorem and de- 
pend only on the isolating integrals of the motion (see e.g.. 



Binney & Tremaine 1987, p. 220). 

Let us assume that the candidate microlensing events 
are caused by Machos in the halo of the Milky Way. This 
remains the most plausible location, although a number of 
investigators (Could, Miralda-Escude & Bahcall 1993; Giu- 
dice, Mollerach, & Roulet 1994) have recently suggested that 
thin or thick discs and spheroids of dark matter objects are 
viable alternatives. Given our hypothesis, the most impor- 
tant thing is to reproduce the structure of the dark matter 
halo as accurately as possible at the locations where typi- 
cal microlensing events occur - between 10 — 40 kpc along 
the line-of-sight from the observer to the Large Magellanic 
Cloud {i = 280°, fe = -33°). Our model of the Milky Way 
has just two components - a luminous thin exponential disc 
and a fiattened dark matter halo. This follows the lead pro- 
vided by Gilmore, Wyse & Kuijken (1989), who claim that 
there is no evidence for dark matter in the disc. The density 
of the bulge and spheroidal component is insignificant at the 
radii of interest and can be neglected. But, it is important 
to include the disc - its contribution to the rotation curve at 
the solar radius is ~ 40% and this alters the normalisation 
of the halo! 

Let us adopt the simple and reasonable assumption that 
the galactic rotation curve is fairly fiat out to the Large Mag- 
ellanic Cloud and deconvolve this into contributions from 
the visible disc and dark matter halo. From stellar kine- 
matic analyses, the local surface density of observable disc 
matter is ~ 48 Mq pc~^ (Gilmore, Wyse & Kuijken 1989). 
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The Milky Way disc - like those of external spiral galaxies 
- is exponential in radius. The scale-length is uncertain, 
but a typical estimate is 3.5 kpc (de Vaucouleurs & Pence, 
1978). 

The equipotentials of dark matter haloes are roughly 
stratified on similar concentric spheroids. A good approxi- 
mation to their gravitational potential is given by (Evans, 
1994) 
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where [R, z) are cylindrical polar coordinates. Here, Re is 
the core radius of the halo and q is the axis ratio of the 
equipotentials. So, q = 1 corresponds to a spherical halo, 
while q < 1 implies the mass model is flattened. The pa- 
rameter vq is a velocity that measures the central depth of 
the potential well. The halo mass density is 
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which falls off like r ^ at large radii. The models - 
known as the power-law haloes - are useful because the 
self-consistent phase space distribution function is known 
analytically and has the simple form 
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where A, B and C are constants given in Evans (1994). Note 
that, in accord with Jeans' theorem, the distribution func- 
tion depends on the two isolating integrals of motion, the 
binding energy E and the angular momentum component 
parallel to the symmetry axis Lz ■ Strictly speaking, the dis- 
tribution of velocities that builds the halo in the combined 
potential field of both disc and halo is what we need - how- 
ever, the potential of the disc is negligible compared to that 
of the halo at the heights above the galactic plane where 
typical microlensing events occur and so the self-consistent 
solution (2) is an excellent approximation. For example, at 
a typical microlens location of half-way along the line-of- 
sight towards the Large Magellanic Cloud, the percentage 
error in neglecting the contribution of the disc to the total 
gravitational potential is < 4%. 

Although the core radius of the halo Re is not well- 
known, our results are insensitive to this parameter. A rea- 
sonable value is ~ 2 kpc (Bahcall, Schmidt & Soneira 1983). 
The axis ratio of the equipotentials q controls the flattening 
of the model. Evidence on the ellipticity of the dark halo is 
sketchy. All we know for sure is that A'^-body simulations 
of gravitational collapse (e.g., Dubinski & Carlberg 1991) 
invariably produce flattened dark haloes. We shall give re- 
sults for the extremes of EO and E6 haloes. The galactic 
constants - the solar radius Rq and the local circular speed 
Veire{Ro) ~ ^^^ uncertain to at least 10%. Merrifield (1992) 
analysed the data on the thickness of the HI layer and con- 
cluded Rq — 7.9 kpc and Veirc(Ro) ~ 200 kms~ . This is 
in accord with a number of recent investigations (Rohlfs et 
al 1986; Reid 1989) which have reported that the lAU es- 
timates (Kerr & Lynden-Bell 1986) of Rq = 8.5 kpc and 
Veire{Ro) = 220 kms~ are too high. Note that increasing 
the local circular speed implies a larger halo density and 
so a greater number of microlensing events. As part of our 



aim is to estimate the minimum number of detections con- 
sistent with a halo built from Machos, it makes sense to 
adopt Merrifield's data-points. The remaining parameters 
Vq and /3 in the power-law halo (1) are chosen to ensure that 
the combined rotation curve of disc and halo is flattish and 
~ 200 kms^ outwards from the solar circle. This yields 
Vq = 138 kms"^ and /3 = -0.2. 

So, our Milky Way model has a luminous thin exponen- 
tial disc embedded in a flattened dark matter halo with a 
simple distribution function. The mass density of the dark 
halo varies like distance" at large radii (c.f., Bahcall, 
Schmidt & Soneira 1983). The model balances simplicity 
with realism. 

3 THE MICROLENSING RATES AND THE 
FREQUENCY DISTRIBUTIONS 

Our aim in this section is to evaluate the observables - the 
microlensing rate F and average duration {te} - together 
with the frequency distributions of events. In the light of 
these calculations, we analyse the data on the possible de- 
tections in Alcock et al (1993) and Aubourg et al (1993). 

The amplification A of the stellar image is related to the 
impact parameter u of the Macho in units of the Einstein 
ring radius R^ by (e.g., Refsdal 1964; Paczyiiski 1986) 
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Let us assume that an observing program can detect events 
above a threshold amplification Ax and impact parameter 
ux- Then, an event occurs whenever a Macho lies within 
a microlensing tube with circular cross-section of radius 
utRe that connects the eye of the observer to the source 
of radiation in the Large Magellanic Cloud. The fatness of 
the tube varies because the threshold impact parameter for 
the deflecting Macho to cause detectable image amplifica- 
tion changes with distance along the line-of-sight. The rate 
is the number of Machos per unit time entering the tube. 
In reality, the tube sweeps slowly through the Milky Way 
as both the observer and the source are in motion. Griest 
(1991) examined the effects of this relative motion and con- 
cluded that it leads to only a slight increment in the total 
rate for the microlensing of Large Magellanic Cloud stars. 
This justifies the approximation of stationary observer and 
source that is used below. 

Let us suppose the radiation source is at a distance L 
from the sun and has galactic coordinates {£, b). By consid- 
ering an element of the microlensing tube of angular extent 
uxREdct and length ds along the line-of-sight, the total 
rate is (Griest 1991) 
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where M is the typical mass of the Machos and VrdvrdvgdO 
is the volume element in cylindrical polar coordinates in 
velocity space. For the Large Magellanic Cloud, the rate 
is 1.35 X 1O"^ut/(M/M0)1/2 event yr'^ if the Milky Way 
halo is spherical, and 1.41x 1O~^mt/(M/M0)1/2 event yx~^ 
if it is as flat as E6. The second observable is the average 
duration of events (te). This is the ratio of the number of 
Machos in the microlensing tube r to the rate at which they 
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enter F, i.e.. 



(ie) = 



(6) 



Now, the number of Machos in the tube - or the optical 
depth - is 



— / nuTREPds. 



(7) 



Again for the Large Magellanic Cloud, the average event 
duration is OMSuriM/MQ)'^/'^ yr for an EO halo, and 
0.368ut{M/Mq)^''^ yr for an E6 halo. Exact expressions 
for the microlensing rate and optical depth in the power- 
law haloes are given in Griest et al (in preparation). 

Let us emphasise that it is not possible to deduce un- 
ambiguously the typical mass of a deflector from the data 
presently available. The best that can be done is to esti- 
mate the most probable mass M that caused an event of 
duration te above the threshold ux- The frequency distri- 
bution of events PMi^e) ~ dT/dte is evaluated by changing 
integration variables in (5) from 6 to te 
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The normalised distribution of events as a function of dura- 
tion is shown in figure 2. The graph is plotted for a threshold 
urp — 1 and Macho mass of O.OSMq. For other values, the te 
axis is scaled by u^M ' and the rate axis by l/iu^M ' ). 
As Griest (1991) argued, the most likely mass giving rise to 
an event of typical duration (te) is the value of M for which 
PAli'te) is largest. This enables us to generate a curve of rel- 
ative probabilities of Macho masses that might have caused 
an event of prescribed duration. 

The candidate event detected by Alcock ot al (1993) 
had a duration of te ~ 0.092 yr and a peak amplification of 
A = 6.86. The high amplification implies that the Macho 
passed close to the centre of the microlensing tube - the im- 
pact parameter is 0.15Re- If w denotes the Macho velocity, 
then the duration of the event is 1.98Re/v- The typical du- 
ration {te) is an average over all amplifications greater than 
the threshold and so is ~ 7ri?£;/(2«). Therefore, the aver- 
age duration (te) implied by the candidate event is roughly 
0.073 yr. (This effect was first noted in Monte Carlo sim- 
ulations by W. Sutherland (1993, private communication)). 
Figure 3 shows the relative probability of Macho masses for 
a spherical halo (full line) and an E6 halo (dashed line), 
from which we deduce the most probable mass of ~ 0.04 or 
0.03 Mq depending on whether the halo is EO or E6. Masses 
of 0.01 Mq and 0.15 Mq are roughly half as likely. 

Now, the end of the main sequence is generally taken 
as 0.08 Mq (e.g.. Nelson, 1990). Objects with masses larger 
than this are dim stars, objects smaller are brown dwarfs. 
The fates of these two classes of object are quite difi'erent - 
low mass stars achieve a state of almost complete thermal 
equilibrium on the order of a Hubble time, whereas brown 
dwarfs are destined to cool to a fully degenerate configu- 
ration. The only fair conclusion from figure 3 is that the 
cause of the candidate event reported by Alcock et al (1993) 
might be microlensing by a dark halo object with a mass in 
the range 0.01—0.15 Mq. This is consistent with both brown 
dwarfs and low mass stars. It is substantially smaller than 
the 0.03 - 0.5 Mq range reported by Alcock et al (1993). 



The EROS collaboration uses a different definition of 
characteristic time. To analyse their data, we first work 
out the timescale over which the amplification exceeded 
A — 1.34 for their events. This gives durations of 0.134 yr 
and 0.154 yr respectively for their first and second detec- 
tion possibilities. Again, we now use the reported amplifica- 
tions of 2.51 and 3.02 to convert to average durations {te) of 
0.116 yr and 0.129 yr respectively. The first event is consis- 
tent with a range of 0.025 — 0.35 Mq, with the most likely 
value ~ 0.10 Mq. For the second event, the most likely 
mass is ~ 0.12 Mq, with masses of 0.03 Mq and 0.45 M© 
half as likely. 

Our final aim is to estimate the minimum number of 
detections that the microlensing programs should find - as- 
suming the halo is completely composed of Machos. An im- 
portant source of uncertainty lies in the value of the local cir- 
cular speed Vcirc{Ro)y which alters the normalisation of the 
halo model. If a higher circular speed is postulated, then the 
density of Machos required to maintain this value out to the 
Large Magellanic Cloud is correspondingly greater, implying 
an increased number of detections. This is why we have cho- 
sen the estimate of Merrifield (1992), which is amongst the 
lowest of recent determinations. Let us assume that the Ma- 
chos have typical masses of the order mass ~ 0.08 Mq . The 
team of Alcock et al (1991) reported monitoring 1.8 x 10 
stars, and so the annual number of detections in the limit 
of 100% efficiency is ~ 8 or 9, depending on whether the 
halo is round or E6. This is the smallest conceivable rate 
consistent with solely baryonic dark matter. This conclu- 
sion depends strongly on our assumption of the absence of 
dark matter in the disc, and more weakly on the errors in 
the Galactic constants (see Griest at al (in preparation) for 
a longer discussion of the microlensing rate scatter) . 

4 CONCLUSIONS 

Microlensing observables are determined by the structure 
of the halo in velocity space. So, we have provided a sim- 
ple and accurate disc-halo model of the Milky Way - with 
distribution functions for the dark matter. This is a good 
representation of the known structure of our Galaxy. How- 
ever, we caution that other models - with rather different 
lensing properties, such as thick discs of dark matter - are 
also compatible with our existing knowledge. 

Our inferred masses for the lensing objects are less than 
those reported by other investigators (see e.g., Alcock et al 
1993, Aubourg et al 1993 and Jetzner & Masso 1993). This 
has a simple explanation. In our model, part of the local 
centrifugal balance is provided by the gravity field of the 
disc. The potential deep in the halo is lowered in compari- 
son with models like the isothermal sphere, where the local 
centrifugal balance is provided by the halo alone. A lower 
potential means that the Machos move on average slower, 
and so smaller masses are needed to give the observed event 
durations. Our most likely values for the masses straddle the 
boundary between low mass stars and brown dwarfs. The 
data is consistent with microlensing by both possibilities. 
Stronger conclusions nmst await more detections. 

If the halo is completely composed of objects of mass 
~ 0.08 Mq, the expected annual number of detections for a 
progam monitoring 1.8 x 10 stars is ^^^ 8 — 9 in the limit of 
100% efficiency. The actual efficiency of the detectors is un- 
der investigation at the moment by the MACHO group. The 
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EROS collaboration have estimated the efficiency of their 
equipment to be approximately 50% (Aubourg et al 1993). 
This suggests that the annual rate of detections should be 
at least 4 or so if the halo is entirely built of Machos. 
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5 FIGURES 

Figure 1: Plot of the contributions of the 
exponential disc and the power— law halo to the 
total rotation curve of our model of the Milky Way 
Galaxy 

Figure 2: The normalised distribution of 
microlensing event rate as a function of event 
duration for a source in the Large Magellanic 
Cloud. The deflecting lens is taken to be a brown 
dwarf or low^ mass star of O.O8M0. The full line 
refers to an EO halo, the broken line to an E6 halo. 
Figure 3: The relative probabilities of Macho 
mass giving rise to events of average duration 
(te) = 0.073 yr. The full line refers to an EO halo, 
the broken line to an E6 halo. 



